home *** CD-ROM | disk | FTP | other *** search
/ Amiga Plus 2000 #5 / Amiga Plus CD - 2000 - No. 5.iso / Tools / Dev / lame_src / mpglib / decode_i386.c < prev    next >
C/C++ Source or Header  |  2000-01-01  |  4KB  |  161 lines
  1. #ifdef HAVEMPGLIB
  2. /*
  3.  * Mpeg Layer-1,2,3 audio decoder
  4.  * ------------------------------
  5.  * copyright (c) 1995,1996,1997 by Michael Hipp, All rights reserved.
  6.  * See also 'README'
  7.  *
  8.  * slighlty optimized for machines without autoincrement/decrement.
  9.  * The performance is highly compiler dependend. Maybe
  10.  * the decode.c version for 'normal' processor may be faster
  11.  * even for Intel processors.
  12.  */
  13.  
  14. #include <stdlib.h>
  15. #if defined(__riscos__) && defined(FPA10)
  16. #include    "ymath.h"
  17. #else
  18. #include    <math.h>
  19. #endif
  20. #include <string.h>
  21.  
  22. #include "mpg123.h"
  23. #include "mpglib.h"
  24.  
  25. extern struct mpstr *gmp;
  26.  
  27.  /* old WRITE_SAMPLE */
  28. #define WRITE_SAMPLE(samples,sum,clip) \
  29.   if( (sum) > 32767.0) { *(samples) = 0x7fff; (clip)++; } \
  30.   else if( (sum) < -32768.0) { *(samples) = -0x8000; (clip)++; } \
  31.   else { *(samples) = ((sum)>0 ? (sum)+0.5 : (sum)-0.5) ; }
  32.  
  33. int synth_1to1_mono(real *bandPtr,unsigned char *samples,int *pnt)
  34. {
  35.   short samples_tmp[64];
  36.   short *tmp1 = samples_tmp;
  37.   int i,ret;
  38.   int pnt1 = 0;
  39.  
  40.   ret = synth_1to1(bandPtr,0,(unsigned char *) samples_tmp,&pnt1);
  41.   samples += *pnt;
  42.  
  43.   for(i=0;i<32;i++) {
  44.     *( (short *) samples) = *tmp1;
  45.     samples += 2;
  46.     tmp1 += 2;
  47.   }
  48.   *pnt += 64;
  49.  
  50.   return ret;
  51. }
  52.  
  53.  
  54. int synth_1to1(real *bandPtr,int channel,unsigned char *out,int *pnt)
  55. {
  56.   static const int step = 2;
  57.   int bo;
  58.   short *samples = (short *) (out + *pnt);
  59.  
  60.   real *b0,(*buf)[0x110];
  61.   int clip = 0;
  62.   int bo1;
  63.  
  64.   bo = gmp->synth_bo;
  65.  
  66.   if(!channel) {
  67.     bo--;
  68.     bo &= 0xf;
  69.     buf = gmp->synth_buffs[0];
  70.   }
  71.   else {
  72.     samples++;
  73.     buf = gmp->synth_buffs[1];
  74.   }
  75.  
  76.   if(bo & 0x1) {
  77.     b0 = buf[0];
  78.     bo1 = bo;
  79.     dct64(buf[1]+((bo+1)&0xf),buf[0]+bo,bandPtr);
  80.   }
  81.   else {
  82.     b0 = buf[1];
  83.     bo1 = bo+1;
  84.     dct64(buf[0]+bo,buf[1]+bo+1,bandPtr);
  85.   }
  86.  
  87.   gmp->synth_bo = bo;
  88.  
  89.   {
  90.     register int j;
  91.     real *window = decwin + 16 - bo1;
  92.  
  93.     for (j=16;j;j--,b0+=0x10,window+=0x20,samples+=step)
  94.     {
  95.       real sum;
  96.       sum  = window[0x0] * b0[0x0];
  97.       sum -= window[0x1] * b0[0x1];
  98.       sum += window[0x2] * b0[0x2];
  99.       sum -= window[0x3] * b0[0x3];
  100.       sum += window[0x4] * b0[0x4];
  101.       sum -= window[0x5] * b0[0x5];
  102.       sum += window[0x6] * b0[0x6];
  103.       sum -= window[0x7] * b0[0x7];
  104.       sum += window[0x8] * b0[0x8];
  105.       sum -= window[0x9] * b0[0x9];
  106.       sum += window[0xA] * b0[0xA];
  107.       sum -= window[0xB] * b0[0xB];
  108.       sum += window[0xC] * b0[0xC];
  109.       sum -= window[0xD] * b0[0xD];
  110.       sum += window[0xE] * b0[0xE];
  111.       sum -= window[0xF] * b0[0xF];
  112.  
  113.       WRITE_SAMPLE(samples,sum,clip);
  114.     }
  115.  
  116.     {
  117.       real sum;
  118.       sum  = window[0x0] * b0[0x0];
  119.       sum += window[0x2] * b0[0x2];
  120.       sum += window[0x4] * b0[0x4];
  121.       sum += window[0x6] * b0[0x6];
  122.       sum += window[0x8] * b0[0x8];
  123.       sum += window[0xA] * b0[0xA];
  124.       sum += window[0xC] * b0[0xC];
  125.       sum += window[0xE] * b0[0xE];
  126.       WRITE_SAMPLE(samples,sum,clip);
  127.       b0-=0x10,window-=0x20,samples+=step;
  128.     }
  129.     window += bo1<<1;
  130.  
  131.     for (j=15;j;j--,b0-=0x10,window-=0x20,samples+=step)
  132.     {
  133.       real sum;
  134.       sum = -window[-0x1] * b0[0x0];
  135.       sum -= window[-0x2] * b0[0x1];
  136.       sum -= window[-0x3] * b0[0x2];
  137.       sum -= window[-0x4] * b0[0x3];
  138.       sum -= window[-0x5] * b0[0x4];
  139.       sum -= window[-0x6] * b0[0x5];
  140.       sum -= window[-0x7] * b0[0x6];
  141.       sum -= window[-0x8] * b0[0x7];
  142.       sum -= window[-0x9] * b0[0x8];
  143.       sum -= window[-0xA] * b0[0x9];
  144.       sum -= window[-0xB] * b0[0xA];
  145.       sum -= window[-0xC] * b0[0xB];
  146.       sum -= window[-0xD] * b0[0xC];
  147.       sum -= window[-0xE] * b0[0xD];
  148.       sum -= window[-0xF] * b0[0xE];
  149.       sum -= window[-0x0] * b0[0xF];
  150.  
  151.       WRITE_SAMPLE(samples,sum,clip);
  152.     }
  153.   }
  154.   *pnt += 128;
  155.  
  156.   return clip;
  157. }
  158.  
  159.  
  160. #endif
  161.